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ABSTRACT 

The occultation of background stellar sources by foreground Kuiper Belt Objects (KBOs) can be used 
to survey physical properties of the KBO population. We discuss statistics related to a KBO occultation 
survey, such as the event duration distribution, and suggest that occultation searches can be effectively 
used to probe the KBO size distribution below 10 km. In particular, we suggest that occultation surveys 
may be best suited to search for a turnover radius in the KBO size distribution due to collisions between 
small-size objects. For occultation surveys that monitor stellar sources near the ecliptic over a few square 
degrees, with time sampling intervals of order 0.1 sec and sensitivity to flux variations of a few percent 
or more, a turnover radius between 0.1 and 1.0 km can be probed. While occultation surveys will probe 
the low-radius limit and imaging surveys will detect KBOs of size 100 km or more, statistics of objects 
with sizes in the intermediate range of around 1 km to 100 km will likely remain unattainable. 

Subject headings: occultations — Kuiper Belt objects — minor planets — solar system: formation 



1. INTRODUCTION 

The presence of minor bodies beyond the orbit of Nep- 
tune is now well established; current estimates suggest that 
there is a total of order 10^ Kuiper Belt Objects (KBOs) 
with sizes greater than 100 km (see Luu & Jewitt 2002 for 
a recent review). Extending KBO population statistics 
to lower sizes than currently known is essential to under- 
standing the formation and evolution of these minor plan- 
ets. Since a kilometer-sized KBO is expected to have a 
magnitude R ^ 30, it is unlikely that imaging surveys will 
be able to obtain reliable statistics of small-size objects in 
the Kuiper Belt. 

The occultation of background stellar sources by fore- 
ground objects has been well utilized as a probe of fore- 
ground source properties (see Elliot & Olkin 1996 for a 
review). In fact, suggestions have already been made to 
detect KBOs via occultation observations (Bailey 1976; 
Brown & Webster 1997; Roques & Moncuquet 2000). On 
the observational side, the Taiwanese-American Occulta- 
tion Survey (TAOS; Liang et al. 2002) plans to conduct an 
occultation search using an array of three to four robotic 
telescopes simultaneously monitoring the same 3 square 
degrees of sky. 

In this Letter, we discuss statistics related to KBO oc- 
cultations and physical properties of the KBO popula- 
tion that can be probed with a wide-field survey targeting 
the ecliptic area. For reasonable parameters that can be 
achieved with ground-based telescopes, we suggest that 
occultation observations are effective at finding small-size 
KBOs with radii between ~ 0.1 and 10 km. At this low 
end of the KBO size distribution, one no longer expects the 
power-law size distribution observed for KBOs larger than 
100 km to hold. The frequent collisions between small-size 
objects, especially at early times during the formation of 
the large bodies, lead to their destruction to dust such that 
a turnover radius is expected at some characteristic size 
scale (Kenyon 2002). For typical values corresponding to 
the Kuiper Belt, with a surface density, E ~ 10^** g cm^^, 
an angular velocity, fl ^ 27r/250 year^^, a KBO density, 
p ~ 1 g cm"-^, and an age for the belt, r 4 x 10^ years. 



an estimate for the turnover radius is Rc ^ tY,{1/ p 0.1 
km. We suggest that occultation surveys have the po- 
tential to either detect or constrain an expected turnover 
radius, which will, in return, help in the theoretical un- 
derstanding of the formation and evolution of the Kuiper 
Belt. 

The discussion is organized as follows. In the next Sec- 
tion, we will summarize important ingredients necessary 
for a calculation of statistics related to KBO occultations. 
In addition to the size distribution of KBOs, we improve 
previous estimates on the KBO occultation rate by pay- 
ing particular attention to the size distribution of back- 
ground stars projected at KBO distances. As we will dis- 
cuss later, the occultation duration distribution is deter- 
mined by both projected star sizes and KBO sizes. In § 3, 
we will present our results and conclude with a summary. 

2. CALCULATION AL INGREDIENTS 

One can characterize an occultation of a background 
star by a foreground KBO with two parameters: the dura- 
tion AT of the occultation and the fractional flux decrease 
A of the background star during the occultation. In terms 
of physical parameters, these two observables are given by 
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respectively. Here, v is the relative velocity of the KBO 
with respect to the observer while R and i?* are respec- 
tively the radii of the KBO object and the background 
star; the latter is the projected stellar radius at the dis- 
tance of the KBO. We have assumed no flux contribution 
from the KBO. In writing the duration of the event, we 
have taken into account the finite size of stellar sources be- 
yond the point-source approximation (cf. Brown & Web- 
ster 1997); we shall see that this is important. 

2.1. KBO distribution 

We make use of the form for the differential size distri- 
bution of KBOs (suggested by Kenyon 2002) as a function 
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Fig. 1. — The distribution of stellar radii projected at a distance 
of 40 AU. We show the distribution as a function of magnitude limit 
in the V-band as labeled on the figure. The left set of four curves 
corresponds to main sequence stars while the right set is for giants. 
The counts are normalized to a total of ~ 56,000 stars (down to 
y ~ 17 — 18) in the three square degree "crowded field" centered 
on I = 170°, 6 = -9° that is planned to be monitored by the TAOS 
survey. 



of radius R such that 
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where Rc is the turnover radius, which is expected to be 
in the range 0.1 to 10 km. The power-law slope, qi, has 
been estimated to be ^ 4 when R > 100 km (Trujillo 
et al. 2001) and we will set gi = 4 hereafter. For the 
present discussion, we treat q2 and Rc as parameters that 
can potentially be probed with occultation observations. 
The normalization constant Aq is set by requiring that 
the total nmnbcr of KBOs above 100 km agrees with the 
estimated number of classical and scattered KBOs, 7 x 10'', 
based on observations so far (Trnjillo et al. 2001). 

The differential distribution of KBOs with distance D. 
df / dD, is taken to scale as with minimum and max- 

imum distances 30 AU and 50 AU respectively (Brown & 
Webster 1997; Kenyon 2002). We set 7 = 2; our results 
are not strongly sensitive to the choice of 7. The KBOs 
are taken to be uniformly distributed in azimuth, with a 
Gaussian distribution P{i) in inclination angle i from the 
ecliptic. Following measurements by Trujillo et al. (2001), 
we take the standard deviation of this Gaussian to be 15°. 

2.2. Stellar distribution 

Following standard approaches discussed in the litera- 
ture, we construct a representative stellar angular radius 
distribution by first creating a population of single stars 
born at a constant rate over the past 10 Gyr, according 
to the initial mass function of Kroupa, Tout & Gilmore 
(1993). This was then evolved to the present day using 
the BSE stellar evolution code (Hurley, Pols & Tout 2001), 
which provides absolute magnitudes and stellar radii as a 
function of time. Finally, the stars were distributed ac- 
cording to a simple exponential Galactic disk model, with 



10=" 



5 



10' 



10" 




10 



10- 10" 
AT (sec) 



10' 




Fig. 2. — The KBO occultation duration (top) and flux change 
(bottom) distributions, with parameters for the KBO size distri- 
bution as shown, for a one-year survey of the same three square 
degrees of sky. The linestyles denote magnitude limits as in Fig. 1. 
In the top plot, we show the variation in the duration distribution 
as a function of 52, the slope of the KBO radius distribution be- 
low the turnover radius, assuming Rc = 1 km. For comparison, we 
also show the duration distribution for a year long survey that tar- 
gets quadrature; while quadrature gives an increase in duration, this 
comes at the expense of a reduction in the total number of events 
by the same factor. In the bottom plot, we show variations to the 
flux change distribution as a function of Rc assuming 92 = 3. 



disk scale radius 3.5 kpc and scale height 200 pc. In- 
corporating extinction corrections from Bahcall & Soneira 
(1980), we selected magnitude-limited samples to produce 
differential stellar angular size distributions, dN/d\ogRi,, 
for a given field, where is the projected stellar radius 
at the distance of the Kuiper Belt. 

The calculated angular size distribution (Fig. 1) rises 
towards small projected sizes, up to a cutoff related to 
the maximum distance to which stars are seen in a given 
sample. The position of the cutoff hence shifts to smaller 
sizes as the magnitude limit is increased. Beyond V ~ 12, 
where the finite extent of the Galaxy becomes important, 
decreasing the Galactic latitude |6| of the field has the same 
effect. The peak at small angular sizes is composed of com- 
parable numbers of A, F, and G main sequence stars. The 
distribution can be observationally obtained for a given 
field, based on colors and magnitudes (Van Belle 1999) 
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and is not likely to be major source of uncertainty. We 
have however neglected the presence of unresolved binary 
companions to half of these stars (Duquennoy & Mayor 
1991). Due to the neglect of light from unresolved com- 
panions, we may have overestimated the number of events 
above a given amplification limit. The neglect of the ex- 
tra cross section presented by companions, however, un- 
derestimates the event rate. Since in most binaries one 
expects the primary to be considerably brighter than the 
secondary, it is likely that the resulting effect is small. We 
expect our predictions for the stellar size distribution to 
be accurate to within a factor of 2. 

2.3. Relative Velocity 

Relative to an Earth-based observer, the velocity of a 
KBO at distance D and observational angle 



is 
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Here, va 



30 km sec and D 



1 AU. As 



is varied, vkbo{D,4>) ranges from at quadrature to 
30(1 — ^JD^/D) km sec~^ at opposition. For a KBO 
at a distance of 40 AU, the maximum velocity is ^ 25 
km sec^^. Since we will concentrate on KBO occultation 
surveys that target areas around the ecliptic, where the 
KBO density is highest, we do not include the dependence 
of KBO relative velocity on inclination angle. For sur- 
veys that monitor specific fields over long periods, such as 
TAOS, note that cj) changes during the year. An alterna- 
tive possibility is to monitor the night sky while keeping 
the position angle constant. 

2.4. Diffraction 

An additional consideration with respect to KBO oc- 
cultations involves diffraction; the Fresnel scale, = 
^yXD/2n, at a distance of 40 AU, for A = 500 nm, IS 
~ 0.7 km and is typical of both KBO and projected stel- 
lar sizes. The resulting effect due to diffraction can be 
described as an increase in the duration of the occulta- 
tion with a corresponding decrease in the flux decrement. 
We therefore modify Eq. (1), by replacing (i? + i?*) with 
-Rdiff, where i?diff is the effective occultation radius in the 
presence of diffraction fringes. We model this following 
Roques & Moncuquet (2002), and assuming that KBOs 
are regularly shaped objects. Note that finite background 
star sizes, especially when i?^ > i?F, reduce effects due 
to diffraction such that a geometrical description is valid 
with i?diff R + R* for long duration events, at the level 
of 0.1 seconds or more. We will highlight the gain involved 
with diffraction, especially for surveys that can reach sen- 
sitivities to flux variations at the level of 1% or below. 

3. OCCULTATION STATISTICS 

With the distributions of R, i?*, and v established 
above, we can now write the expected number of KBO 
occupations during a total observing time Ttot as 
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The integrals are taken such that one allows for a min- 
imum fractional dimming, ^mim a-nd a minimum sam- 
pling interval in time, ATmin, with which flux from stars 
down to some magnitude limit can be monitored. In 
a real survey, Amin will be a function of magnitude, 
and we consider an example of this later in the Section. 
Instead of discussing the total number of occultations, 
we consider the fractional distribution of occultations in 
terms of observables, mainly event duration, fractional 
flux variation, and KBO radius. These are calculated 
by introducing the necessary Dirac delta function before 
taking integrals in Eq. 4, such that, say, dNtot/dST = 
]d<i>...JdR dNKBO /dR2RdisvTtotSD {ST - 2i?diff/u). 

In Fig. 2, we show the expected number of KBO occul- 
tations as a function of duration, for a year-long survey 
of a given field, using our geometrical description for oc- 
cultations, with no corrections due to diffraction. Most of 
the occultation events are of short duration, and are asso- 
ciated with the main sequence stars producing the peak at 
small angular size. For comparison, we also show the du- 
ration distribution towards quadrature. While durations 
increase by an order of magnitude or more, there is an 
equivalent reduction in the total number of events, which 
can be understood by noting that the duration scales as 
while the total number of occultations scales as v. 
For a long term survey of a given field, the distribution 
predicted towards opposition alone is similar to that ob- 
tained when one integrates over position angle throughout 
the year. This is due to the fact that, away from quadra- 
ture, which is highly localized, all other position angles 
lead to comparable velocities. 

The predictions in Fig. 2 reveal that for surveys which 
are limited to large timing intervals, of order 0.1 sec or 
more, but fainter magnitudes, V ~ 16, observations toward 
quadrature are desirable; this, however, requires monitor- 
ing of stars in a number of different fields, since the back- 
ground sky shifts with respect to the direction of quadra- 
UiTC during a long period of observation. With small im- 
provements in the minimum duration to which a survey is 
sensitive, one can obtain a significant increase in event rate 
by observing, say, one field throughout the year, instead of 
targeting the low-velocity quadrature region. The use of 
a single field also allows better control of systematics and 
a reduction of confusions, such as those due to intrinsic 
stellar variabilities. 

The bottom plot of Fig. 2 shows the expected distribu- 
tion of flux variations. This distribution has a power-law 
behavior when A < 0.1 and flattens out as A approaches 
one. Note that, in a strict sense, the maximum value of 
A is one, when a KBO fully occults the background star. 
Here, for illustration purposes, we have allowed A to vary 
beyond one in order to show the flattening of this distri- 
bution. 

In Fig. 2, we also show variations in the expected num- 
ber of observed occultations induced by variations in the 
ciitoff radius Rc, and the slope q2 below the cutoff radius. 
With sufficient sensitivity to flux variations and with a 
high sampling rate, one can also probe these two param- 
eters; for example, the expected number of occultations 
with flux variations between 1% and 10% differs by two 
orders of magnitude when Rc is varied between 0.1 km 
and 1 km. 

In Fig. 3, we show the expected number of occultations 




Fig. 3. — The KBO occultation distribution as a function of the 
KBO size. For occultation surveys with sensitivity to flux variations 
of order 1% and sampUng intervals of order 100 msec, KBOs down to 
a radius of ~ 0.1 km can be probed. For such surveys, one does not 
gain significantly by monitoring stars with magnitudes fainter than 
■~ 14. In the top plot, we also show the expected boost in number 
of events due to diffraction for the best case scenario involving a 
sensitivity to flux variations at a level of 1%. In the bottom plot, 
we consider a single experiment for which FA^ . AT[uin=constant 
with normalization set at = 10% at 14th magnitude with 

^Tlnin = 0.2 sec. We improve ATmjn as shown on the plot, for 
a survey with V < 14, except for the ease of AT^i^ = 0.02 see, 
for which we show curves for limiting magnitudes of 12, 14, and 
16. Few additional events are gained by monitoring deeper than 
V = 14, especially when one considers the additional data reduction 
required. 



as a function of KBO radius. For surveys with sampling 
intervals of 0.1 seconds and sensitivity to normalized flux 
decreases of 1% or better, one probes the size distribu- 
tion down to a radius of 0.1 km. The cutoff corresponds 
to the KBO radius producing the smallest detectable flux 
decrease A^in for a star whose size leads to an occultation 
with AT = ATmin- We also highlight the gain associated 
with diffraction here. With a detection threshold of order 
1% in flux, and a time sampling interval of 0.01 sec. diffrac- 
tion boosts the number of expected occultations by an or- 
der of magnitude. It is unlikely that occultation surveys 
with ground based small telescopes can achieve the high 



flux sensitivity to see effects related to diffraction, though 
targeted surveys of smaller stars with 8m-class telescopes 
may be used to exploit the gain associated with diffraction 
(Roques & Moncuquet 2000). 

We now consider a real experiment: the planned TAOS 
survey (Liang et al. 2002). We take into account the pres- 
ence of photon noise, which leads to a decrease in flux 
sensitivity when individual exposures are shorter. There 
is a tradeoff between sensitivity and sampling rate for 
a given experimental setup, which can be expressed as 
A^jjjFATmin = constant, where F is stellar flux. In Fig. 3 
bottom plot, we show the expected distribution of occul- 
tations as a function of the KBO radius. Here, we as- 
sume our fiducial KBO population {qi = 4, (72 = 3 and 
Rc = 0.1 km) and normalize the scaling relation with the 
expected TAOS values of ATmin = 0.2 sec, and ^min — 
10% at 14th magnitude. As shown in the bottom plot of 
Fig. 3, with improvements to the sampling interval, and 
a corresponding decrease in the flux sensitivity, one low- 
ers the minimum KBO radius probed. Following Fig 2, 
this is due to the rapid increase in the number of events 
with decreasing duration, when compared with the distri- 
bution in terms of flux variations, which is slowly varying 
when A ~ 10%. We see that for stirvcys like TAOS, sen- 
sitive to flux variations of order 10% at a magnitude of 
14 and time intervals of a few tenths of a second, there is 
no significant gain in monitoring stars fainter than 14th 
magnitude, since a large fraction of the detectable occul- 
tations are associated with bright stars of large angular 
size, which have smaller mimimum detectable flux varia- 
tions and larger occultation durations, so more detectable 
occultations. The additional data procesing involved with 
monitoring fainter stars could also be cumbersome, given 
the diminishing returns. 

To summarize, occultation surveys are ideal for prob- 
ing a turnover radius in the size distribution expected in 
certain theoretical models of formation and evolution of 
KBOs. To obtain adequate statistics, ground-based sur- 
veys with smaller telescopes should monitor bright stars 
{V ^ 14) on and around the ecliptic, in fields where the 
density of stars is also high. In the case of the planned 
TAOS survey (Liang et al. 2002), it may be advisable to 
optimize; the survey by decreasing the minimum time in- 
terval with a corresponding decrease in sensitivity to fiux 
variations. Beyond the TAOS survey, we encourage fur- 
ther occultation surveys to probe the low end of the size 
distribution since this is where theoretical models of KBOs 
can be tested in detail. 
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